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OPTICAL AIR DATA SYSTEMS AND 
METHODS 

RELATED APPLICATIONS 

5 

This application is a continuation and claims benefit of 
priority to commonly -owned and copending U.S. patent 
application Ser. No. 1 1/103,020 filed Apr. 1 1, 2005 which is 
a continuation ofU.S. patent application Ser. No. 10/632,735, 
filedAug. 1, 2003, now U.S. Pat. No. 6,894,768, which claims to 
benefit of priority to U.S. Provisional Patent Application No. 
60/400,462, filedAug. 2, 2002. All of the above-mentioned 
applications are incorporated herein by reference. 

U.S. GOVERNMENT RIGHTS 15 

This invention was made in part with the support of the 
U.S. Government; the U.S. Government has certain rights in 
this invention as provided for by the terms of Grant # NAS4- 
0243 awarded by the NASA Dry den Flight Research Center. 20 

BACKGROUND 

An Air Data System (‘ADS”) provides sensed telemetry 
informing pilots, navigators or Vehicle Management System 25 
computers of air parameter(s) affecting aircraft stability. 
These air parameters include, for example, air speed, air 
temperature and air pressure, each being useful for navigation 
and flight control. The ADS exists in many forms, for 
example, as mechanical, opto -mechanical or opto-electronic 30 
devices. 

One mechanical ADS includes a Pitot tube, which affixes to 
the exterior of the aircraft. The Pitot tube is a pneumatic 
measuring instrument that receives air external to the aircraft 
and determines air speed based upon air pressure. A typical 35 
Pitot tube uses a pressure transducer to measure the pressure 
of air received through a plurality of holes. Determined air 
pressure is combined with air density measurements and air 
temperature measurements using the well-known Bernoulli 
equation to calculate air speed. Pilots and Vehicle Manage- 40 
ment Systems use this information for decisions regarding 
aircraft stability, flight control and navigation. 

The prior art ADS has limitations. With respect to the Pitot 
tube, for example, if aircraft velocity is too low, the pressure 
transducer may lack requisite sensitivity to provide meaning- 45 
ful telemetry data, such that transducer signal noise or error 
may be greater than meaningful differences in pressure mea- 
surement. Alternatively, if the velocity is very high (e.g., 
supersonic), certain assumptions, such as those regarding 
incompressibility of air flow in Bernoulli’ s equation, are vio- 50 
lated and the measurement is incorrect. Accordingly, Ber- 
noulli’s equation must be altered to compensate for various 
breached assumptions. This compensation process must con- 
tinue in rapidly changing conditions of altitude and air den- 
sity. It is, therefore, desirable to improve accuracy and oper- 55 
ability of the ADS. 

One improvement is an Optical Air Data System 
(“OADS”), which uses light to determine parameters of air 
speed. The OADS transmits light pulses into the atmosphere 
and receives light that is reflected, or “backscattered”, from 60 
aerosols towards the aircraft. Aerosols are fine solids and/or 
liquid particles suspended in air or other gases. 

The OADS may also measure the Doppler effect, receiving 
the backscattered light and measuring return frequency to 
determine speed. However, the prior art OADS relies on 65 
scattered light that is unpredictable because of varying aero- 
sol distributions. For example, aerosol distribution varies sig- 
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nificantly with altitude and cloud content. In addition, some 
regions of the atmosphere contain too few aerosols to enable 
reliable air data measurements. Thus, aerosol-based optical 
air data systems cannot determine the air speed at all altitudes 
frequented by modern aircraft. Finally, an OADS that relies 
solely on aerosol scattering cannot determine air temperature 
or air pressure. Air temperature and pressure are critical air 
data parameters for determining air density and Mach num- 
ber. Hence, there is a need to improve optical air data systems 
and methods. 

SUMMARY OF THE INVENTION 

In one embodiment, an Optical Air Data System (“OADS”) 
advances the art and overcomes the problems mentioned 
above by determining air parameters (air speed, air tempera- 
ture and air pressure) based upon laser energy that is back- 
scattered from both air molecules and aerosols. In one 
example, the OADS projects laser energy to the air and 
receives backscattered laser energy from the air molecules 
and/or aerosols. The OADS converts such backscattered laser 
energy to electronic signals which are processed to determine 
desired air parameters. Such air parameters may include, for 
example, air speed, air temperature and air pressure, among 
others, and provide telemetry for air vehicle flight control and 
navigation. 

In one embodiment, the OADS includes at least one laser 
configured for generating laser energy. The OADS also 
includes one or more transceivers configured for projecting 
the laser energy to the air. Each transceiver subsequently 
receives backscattered laser energy from the air. A portion of 
the backscattered laser energy received by each transceiver is 
filtered prior to electronic signal conversion. 

In one embodiment, each transceiver includes a vapor filter 
that filters the backscattered laser energy within a predeter- 
mined band. In one example, a transceiver includes a fixed 
frequency vapor filter, such as an atomic vapor filter or a 
mercury vapor filter, that provides “notch” or “band stop” 
filtering of the backscattered laser energy. The band stop filter 
may provide peak absorption of background light within the 
band stop region of the filter. In one embodiment, the laser is 
tunable and generates the laser energy at a center-tuned wave- 
length that substantially corresponds to a center wavelength 
of the band stop of the vapor filter. In another embodiment, 
the center wavelength of the filter corresponds to an absorp- 
tion region for ultraviolet wavelengths occurring in ozone 
(0 3 ). Once filtered, a transceiver converts the backscattered 
laser energy into a computer proces sable electronic signal for 
determination of the desired air parameters. 

In one embodiment, each transceiver includes one or more 
detectors or transducers configured to receive filtered and 
unfiltered backscattered laser energy for conversion to elec- 
tronic signals. In one example, a transceiver includes an opti- 
cal transducer that converts a particular range of wavelengths 
to computer processable signals used to determine the desired 
air parameters. Accordingly, in one embodiment, the OADS 
includes a computer for processing the electronic signals to 
determine the desired air parameters. Other features are 
apparent within the description that follows. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 shows one Optical Air Data System (“OADS”); 

FIG. 2 shows one OADS; 

FIG. 3 shows one graph useful in illustrating exemplary an 
air speed calculation with an OADS; 
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FIGS. 4-9 show graphs illustrating exemplary calculations 
for other air parameters with an OADS; 

FIG. 10 shows one OADS transceiver; and 

FIG. 11 shows a flowchart of one exemplary methodical 
operation of an OADS. 5 

DETAILED DESCRIPTION OF THE DRAWINGS 

FIG. 1 shows one Optical Air Data System (“OADS”) 101 
mounted on or within an aircraft 102. In this embodiment, io 
OADS 101 is configured for projecting laser energy 103 (e.g., 
laser radiation) to air 104. Impingement of laser energy 103 
on air 104 and aerosol particles 105 (in air 104) causes scat- 
tering of laser energy 103 , which is represented in FIG. 1 as an 
energy scatter field 106. Distance between aircraft 102 and 15 
energy scatter field 1 06 is controlled by the amount of overlap 
between the out-going laser energy 103 and the transceiver 
110 field of view at a distance from aircraft 102 , to provide an 
optimized intensity for return laser energy 107 and to elimi- 
nate the overwhelming measurement error arising from dis- 20 
placed air proximate aircraft 102. OADS 101 detects return 
laser energy 107 that is backscattered from air 104 at laser 
energy scatter field 106. Laser energy 107 may be in the 
ultra-violet (UV) spectrum, for example, having a wave- 
length within a range of 250 nm to 270 nm. Other ranges 25 
operable to produce scatter field 1 06 may be used as matter of 
design choice. 

Return laser energy 107 typically contains molecular scat- 
tered (e.g., Rayleigh) components 107A and/or aerosol scat- 
tered (Mie) components 107B. OADS 101 distinguishes the 30 
molecular scattered components 107A from the aerosol scat- 
tered components 107B and correspondingly determines one 
or more air parameters based on the scattered laser energy 
107. Examples of such air parameters include air speed, air 
pressure, air temperature and aircraft orientation angles rela- 35 
five to the local wind. OADS 101 may be configured with 
other aircraft as well, such as unmanned air vehicles (UAVs), 
helicopters, gliders and space shuttles. Although illustrated 
within a “nose” 108 of aircraft 102, OADS 101 may be 
configured in any other areas of aircraft 102 . 40 

As shown in FIG. 1, OADS 101 includes a laser 109 con- 
figured for generating laser energy 103. Transceiver 110 is 
configured for transmitting laser energy 103, from laser 109 
via optical coupling 111 , and receiving backscattered laser 
energy 107. Optical coupling 111 may exist in the form of a 45 
fiber optic connection or free space transmission. Accord- 
ingly, transceiver 110 projects the laser energy as laser radia- 
tion 103 to air 104. Air 104 scatters laser radiation 103 at 
energy scatter field 106 in a plurality of directions (e.g., 
illustrated as vectors 112). Scatter field 106 also returns, or 50 
backscatters, laser radiation 107 towards transceiver 110, 
which subsequently receives the backscattered laser radiation 
107. Transceiver 110 converts backscattered laser radiation 
107 to processable electronic signals, via computer 113, to 
determine the air parameters. 55 

Computer 113 communicatively couples with transceiver 
110 and processes signals from transceiver 110 to distinguish 
molecular scattered component 107A from aerosol scattered 
component 107B. Computer 113 determines the air param- 
eters based on laser radiation 107 backscattered from mol- 60 
ecules and/or aerosols in air 104. Accordingly, as described 
below, computer 113 employs one or more digital signal 
processing algorithms to determine such parameters. 

While OADS 101 illustrates one transceiver 110 in an 
exemplary embodiment, a plurality of transceivers may be 65 
used depending on application. For example, a helicopter 
employing OADS 101 may use two transceivers 110 to deter- 
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mine air parameters such as forward velocity (e.g., air speed) 
and a horizontal plane, or “yaw”, of the helicopter. An air- 
plane may use three transceivers 110 positioned in a particu- 
lar manner to determine various aircraft geometries, such as 
angle of attack and sideslip, in addition to the air parameters 
of air speed, air pressure and air temperature. In addition, air 
vehicles (fixed wing and rotary) may employ three or more 
transceivers and/or lasers to increase Air Data System reli- 
ability through a redundant system architecture. Using three 
OADS transceivers mounted orthogonally to one another 
may fully resolve a total airspeed vector by providing three 
independent measurements for the air speed vector (i.e., cor- 
responding to three axes of a Cartesian coordinate system). 
The transceivers are located in uncommon planes and their 
geometry respective of an aircraft center-line known. Vector 
algebra may then be used to determine the hill airspeed vec- 
tor, including forward air speed, angle-of-sideslip and angle- 
of-attack. 

FIG. 2 shows one OADS 140. OADS 140 illustrates 
another embodiment used for determining air parameters, 
such as those described in FIG. 1, based upon laser energy 
backscattered from both air molecules and aerosols. In this 
embodiment, OADS 140 includes laser 141 configured for 
generating laser energy 142. Laser 141 may be tunable laser 
having a tuned center wavelength of about 253.7 nm, 
although other wavelengths may be used. For example, laser 
141 may be a frequency quadrupled, Nd:YAG (i.e., neody- 
mium :yttrium-aluminum-gamet) pumped Ti: Sapphire (tita- 
nium-sapphire) laser. Alternatively, frequency-quadrupled 
Yb-doped (ytterbium- doped) fiber lasers can be used that 
offer important benefits of smaller size, lighter weight, 
increased robustness and improved reliability, as compared to 
Nd:YAG-pumped Ti: sapphire lasers. Laser 141 generates 
laser energy tunable across 40 GHz in about 100 MHz spec- 
tral steps. Additionally, laser 141 may tune +/-20 GHz about 
a center frequency of approximately 1 182.5 THz, or c/253.7 
nm, where c is the speed of light (approximately 3x1 0 8 m/s). 
Laser 141 radiates laser energy 142 to beam splitter 143, 
which splits the beam into two components 143A/143B. 
Component 143A is directed through air 144; component 
143B is directed to beam splitter 145. 

In particular, component 143 A of laser energy 142 directed 
to air 144 is scattered into scatter field 146. Scattering of 
component 143A is illustrated by scattering vectors 147 in 
scatter field 146, whereas return scattering is illustrated by 
backscattered laser radiation 148. Component 143 B of the 
laser energy 142 is used as a reference for comparison to 
backscattered laser radiation 148. Such a comparison is for 
example useful in determining air parameters such as air 
speed, since transmitted and received frequencies of the laser 
energy can be ascertained for use in a Doppler equation; such 
a process is explained in greater detail herein below. 

In the illustrated embodiment, backscattered laser radia- 
tion 148 is received through optics 149. In one example, 
optics 149 is a telescope that demagnifies backscattered laser 
radiation 148 into a laser beam 150. Optics 149 also directs 
beam 150 to beam splitter 151, to split beam 150 into two 
components 150A/150B. Component 150B of beam 150 
passes through vapor filter 152 to detector 153 to produce 
electronic signal 158 representative of the component 150B 
impinging detector 153 ; whereas component 150A is directed 
by beam splitter 151 to detector 154. 

In one embodiment, detector 154 is a photodiode that 
receives laser 148 radiation and converts such radiation into 
an electronic signal 155. Detector 154 connects to a central 
computer 156 to process electronic signal 155. Similarly, 
detector 153 is a photodiode configured for detecting com- 
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ponent 150B, which is filtered by vapor filter 152 as filtered 
component 157. Detector 153 converts component 157 to an 
electronic signal 158 for processing by central computer 156. 

Accordingly, electronic signal 158 corresponds to back- 
scattered laser energy 148 as filtered by vapor filter 152; and 
electronic signal 155 corresponds to unfiltered backscattered 
laser energy 148. Electronic signal 155 is thus used to nullify 
certain anomalies as computer 156 processes electronic sig- 
nal 158. For example, when processed with electronic signal 
158, signal 155 may be used to remove, from signal 158, 
certain laser transmission power fluctuations in filtered com- 
ponent 157 caused by atmospheric changes in air 144. Such a 
process is explained in more detail in FIGS. 4-7. 

Reference component 143B of the laser energy 142 is split 
into two components 159/160 by beam splitter 145. Compo- 
nent 1 6 0 is directed by beam splitter 1 45 to vapor filter 1 52 via 
mirrored surface 161, to measure filter characteristics, 
whereas component 159 is directed by beam splitter 145 to 
detector 162, to generate electronic signal 163. Electronic 
signal 163 is for example used to normalize power fluctua- 
tions in the return of backscattered laser radiation 148 caused 
by power fluctuations in the generation of laser energy 142 by 
laser 141. Such a process is explained in more detail in FIGS. 
4-7. 

Vapor filter 152 filters component 160 to produce filtered 
component 164. Filtered component 164 is directed to detec- 
tor 165, via mirrored surface 166, and then converted to an 
electronic signal 167. Central computer 156 processes elec- 
tronic signal 167 to determine filter characteristics, such as 
frequencies and suppression features of the band stop region 
of vapor filter 152. Such a process is also explained in more 
detail in FIGS. 4-7. 

It should be noted that FIG. 2 displays the OADS 140 using 
free space optical transmission and optical components, such 
as beam splitters 143, 145 and 151, and mirrors 161 and 166. 
However, a more environmentally robust method may use 
optical fiber for laser 141 energy transmission along paths 
142, 143A, 143B, 159, 160, 164, 150, 150A, 150B and 157. In 
such an embodiment, fiber splitters are used for 143, 151 and 
145 and mirrors 161 and 166 are eliminated. 

FIG. 3 shows one graph 200 useful in illustrating an exem- 
plary air speed calculation with OADS 140. Graph 200 shows 
two curves 201/202 comparing normalized laser energy mag- 
nitudes as a function of frequency (normalized laser energy 
magnitude plotted with respect to axis 205 and frequency 
plotted with respect to axis 204). Curve 202 exemplifies fil- 
tered radiated laser energy such as that of filtered component 
164 of FIG. 2. As such, curve 202 shows filter characteristics 
of vapor filter 152 of FIG. 2 determined by processing of 
electronic signal 167. Curve 202 shows a peak absorption of 
filter 152 occurring at a down- translated frequency of 0 GHz. 
By way of example, the actual peak absorption frequency of 
filter 152 may be about 1182.5 THz (i.e., having a corre- 
sponding wavelength of about 253.7 nm). 

Laser energy 142 generated by laser 141 passes through 
filter 152 to provide filtered component 164. Once filtered 
component 164 is converted to electronic signal 167 by detec- 
tor 165, computer 156 analyzes and stores features of vapor 
filter 152 through digital signal processing of signal 172 (e.g., 
reference features under controlled conditions for use in 
future calculations). As shown in this example, features of 
vapor filter 152 have approximately 10% normalized absorp- 
tion at approximately +1-5 GHz (i.e., 0.9 normalized trans- 
mission factor at approximately +/- 5 GHz according to axis 
205) about the peak absorption frequency. Other types of 
suitable filters may include different absorption/transmission 
features. 
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Curve 201 exemplifies filtered backscattered laser energy 
such as that of filtered component 157 of FIG. 2. In one 
embodiment, curve 201 is used to determine air speed by 
comparison to curve 202. For example, curve 202 illustrates 
5 how vapor filter 152 affects laser energy 142; curve 201 
similarly illustrates how vapor filter 152 affects laser energy 
142 as laser energy 142 is backscattered (e.g., returned laser 
radiation 148) from air 144. Frequency shift 203 represents 
the change in frequency of peak absorption for vapor filter 
to 152 between transmitted laser energy 142 and returned laser 
energy 148. Computer 156 processes algorithms applying 
Doppler velocity equation to determine air speed from fre- 
quency shift 203. 

To determine air speed in one embodiment, computer 156 
15 determines how far in frequency the peak absorption fre- 
quency of filtered component 157 has shifted from the initial 
laser frequency by comparing curve 202 to curve 201 (e.g., 
comparing peak absorption frequencies of filtered compo- 
nents 157/164). The frequency shift 203 between the two 
20 frequencies substantially equates to a radial wind velocity 
through the Doppler velocity equation: 


2 v R (Eq. 1) 


where Av n represents the Doppler frequency shift, \ R rep- 
resents velocity component of the vehicle (e.g., aircraft 101 of 
FIG. 1) along the laser direction of propagation 143 B and A, 
30 represents the wavelength of laser energy 142. 

In one embodiment, the wind velocity can be measured by 
determining the frequency shift from curve 202 of graph 200 
as compared to curve 201 of graph 200. This is accomplished 
by calculating the symmetry point of both spectra 202 and 
35 201 and determining the difference in symmetry points 
between the two curves. 

FIGS. 4-9 show graphs illustrating exemplary calculations 
for other air parameters with OADS 140. For example, after 
determining frequency shift due to air speed as shown in FIG. 
40 3, other air parameters, such as air temperature and air pres- 
sure may be calculated. In one example, computer 156 ini- 
tially determines an intensity measurement of the detected 
backscattered laser radiation (e.g., filtered component 157 
detected by detector 153) from electronic signal 158. This 
45 experimentally verified intensity measurement of returned 
laser energy corresponds to the following equation: 

$s(. v ) = P L?LDsTjiEddvJdv l [L(v l )F(v r -v)(rR(v r - (v 7 - 

Av D ))+wM(v r -Av D )))] (Eq. 2) 

50 where S s (v) is electronic signal 158 from detector 153; P z 
is the laser power, T L is the transmission coefficient through 
air 144 along laser path 143B, L(v ; ) is the laser line shape 
inherent to the laser 141 output as a function of laser fre- 
quency v ; , T^ is the transmission coefficient through air 144 
55 along laser path 148, E^ is optical efficiency of the detector 
channel through detector 153, F(v) is the band stop frequency 
range of vapor filter 152 centered at a frequency of v, R is 
Rayleigh scattering as a function of frequency (applicable to 
the Rayleigh regime) v r for backscattered laser radiation 
60 minus the quantity of laser frequency v ; minus the Doppler 
shift Av^,, r is the Rayleigh scattering magnitude coefficient 
dependent on air density and the Rayleigh backscattering 
coefficient, M is Mie scattering as a function of v r minus the 
quantity of v l minus Av^, m is the Mie scattering magnitude 
65 coefficient dependent on aerosol concentration and the Mie 
backscattering coefficient, and D s is detector 153 efficiency 
of the detector channel 158. The Rayleigh backscattering 



US 7,760,339 B2 


7 

coefficient r and the Mie backscattering coefficient m are 
constant for a particular atmosphere. These coefficients cor- 
respond to the number of scatterers (i.e., molecules for Ray- 
leigh, aerosols for Mie) per unit volume of atmosphere. 

With a result obtained for the measured intensity of the 
returned laser energy, computer 156 commences the process 
of determining the other air parameters. Such a process, for 
example, begins with determining characteristics of vapor 
filter 152 via the transmission of reference laser energy 
through vapor filter 152 (e.g., component 160). For example, 
measuring the band stop characteristics of vapor filter 152 
with laser 141 (e.g., via component 143B to electronic signal 
167) during experimentation yields a convolution of the laser 
wavelength and the filter according to the following equation: 

S^P^DjdvjfLivdFivrv)], (Eq. 3) 

where S^v) is signal 167 from detector 165 as a function of 
frequency v (e.g., as illustrated in curve 122 ); E^is the optical 
efficiency of filter 152 collection along paths 160 and 164, 
and Dp is the optical efficiency of detector 165. 

Note that all optical efficiencies Ep. and E^ capture signal 
losses that are optical in nature. For example, the optical 
efficiency for detector 165, Dp., includes the optical beam 
splitting ratios for beam splitters 143 and 145, the optical 
transmission and coupling across filter 152 and the optical 
delivery efficiency onto detector 165. The optical collection 
efficiency for detector 153, E^, includes the collection effi- 
ciency of telescope 149, the optical coupling efficiency into 
path 150, the beam splitter ratio of beam splitter 151, the 
transmission efficiency across filter 152 and the delivery effi- 
ciency onto detector 153. Detector efficiencies Dp. and D s 
include the detector conversion efficiencies for detectors 165 
and 153, respectively. Thus, Dp. is the conversion efficiency 
whereby detector 165 converts laser energy along path 164 
into an electrical signal 167. Likewise, Ds is the conversion 
efficiency whereby detector 153 converts laser energy along 
path 157 into an electrical signal 158. 

Backscattered laser radiation 148 can include power fluc- 
tuations that are caused by laser 141 while generating laser 
energy 142. Accordingly, laser energy detected by detector 
162 (e.g., via component 159) assists in normalizing power 
fluctuations attributable to laser 141 by removing such power 
fluctuations from the analysis. In one embodiment, detector 
162 converts the laser energy into electronic signal 163. In 
turn, computer 156 processes and normalizes according to the 
following equation: 

S L =P L E L DjdvL(y), (Eq. 4) 

where Sp is the electronic signal 163 from detector 162, E L 
is the optical collection efficiency for detector 162, D z is the 
conversion efficiency of detector 162 and P z is the power of 
laser 141. Note that the optical collection efficiency E L 
includes the beam splitting ratios of beam splitters 143 and 
145 and the delivery efficiency of laser beam path 159 onto 
detector 162. 

Curve 221 of FIG. 4 represents the magnitude of laser 
energy (component 164) filtered by vapor filter 152 and nor- 
malized between 0 and 1 . Curve 221 represents the magnitude 
of the laser energy as a function of frequency (i.e., laser 
energy magnitude plotted with respect to axis 222 and fre- 
quency plotted with respect to on axis 223). Curve 221, there- 
fore, illustrates filtered laser energy via component 160 as 
determined by computer processing of electronic signal 167, 
plotted as laser energy magnitude normalized between 0 and 
1 versus frequency. 

In one embodiment, the absorption/transmission charac- 
teristics of the vapor filter 152 are normalized using Eqs. 3 
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and 4. Eq. 3 yields the stop band characteristics of filter 152 
and Eq. 4 accounts for power fluctuations in the generation of 
laser energy 142. With the power fluctuations of Eq. 4 sub- 
stantially removed, a “normalization channel” is revealed and 
5 power fluctuations attributable to atmospheric changes can be 
accounted for. 

In one embodiment, additional power fluctuations caused 
by atmospheric changes in air 144 are also removed. For 
example, laser energy detected by detector 154 (e.g., via 
to component 150A) assists in removing laser power fluctua- 
tions caused by atmospheric changes in air 144. Accordingly, 
detector 154 converts received laser energy into electronic 
signal 155. Computer 156, in turn, processes electronic signal 
155 to determine the normalized laser energy magnitude 
15 according to following equation: 

Sn= p lT l TrE n Dn$ M cfy/LCv,) (rR (v- (v r A v D ))+mM 

(v-(v r AvD))], (Eq. 5) 

where S^is the signal 155 from detector 154; E A ,is optical 
20 collection efficiency of the detector 154 and D^ is the con- 
version efficiency of detector 154. 

In one embodiment, it is advantageous to normalize the 
various characteristic functions to enable a closed-loop solu- 
tion to the process of determining temperature and pressure. 
25 In one example, therefore, computer 156 calculates the nor- 
malized laser line shape according to following equation: 

jL(vi)dvrl, (Eq. 6) 

where (as before) v I is laser line shape frequency and L 
30 denotes the laser line shape as a function of frequency. In 
another example, computer 156 calculates normalized Ray- 
leigh Function according to the following equation: 

SR(v r )dv= 1, (Eq. 7) 

35 where R denotes the Rayleigh line shape as a function of 
frequency v r , applicable to the Rayleigh regime. In another 
example, computer 156 scales the electronic signal 167 
recorded from detector 1 65 by dividing all recorded values by 
the maximum value according to the following equation: 

40 MAX(5p(v))=l, (Eq. 8) 

where MAX denotes an operation that finds a maximum 
value of a particular function, and denotes electronic sig- 
nal 167 measured from detector 165, as a function of fre- 
quency v (e.g. as illustrated in curve 221 of FIG. 4). In another 
example, computer 156 normalizes the Mie Function accord- 
ing to the following equation: 

M(v)=S(v), (Eq. 9) 

where 6 (v) is the delta function. 

50 • 

In one embodiment, dividing the signal 1 67 collected from 
detector 165 (and represented by Eq. 3, above) by the signal 
163 collected from detector 163 (and represented by Eq. 4, 
above) removes laser 141 power fluctuations, as follows: 

55 

Sp{v) P L E F D F Jd Vl [U Vl )F( Vl - v)] (Eq ' 10) 

SM ~ P L E L D L fdvL(v ) 

60 

Equation 10 can be simplified to: 


Spiy) _ (Eq. 11) 

65 S L (v) E l D l W 
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In one embodiment, tuning the laser 141 to a frequency 
(vj far enough removed from the effects of the vapor filter 
152 enables the measurement of the ratio of the optical and 
detector efficiencies of the signal channels 167 (S^, repre- 
sented by Eq. 3 above) and 163 (S i; represented by Eq. 4 5 
above). This, in turn, enables the normalization of the signal 
167 measurement to one, for simultaneously checking for 
laser, detector and filter abnormalities on a scan-by-scan 
basis: 


Sf(Vq) = E f D f (Eq. 12) 

Sl(v 0 ) ElDl 

In one embodiment, LF(v) are determined to generate a 15 
look up table of the convolution of theoretical Rayleigh func- 
tions (calculated in terms of temperature and pressure) with 
the measured filter function. Since, the measured filter func- 
tion is already the convolution of the laser and filter spectra, ^ 
convolving the Rayleigh function with the measured filter 
signal 167 yields the expected return signal from an atmo- 
sphere of pure Rayleigh scatterers. 

In one embodiment, the measured signal 158, which is the 
backscatter return from the atmosphere 144, that passes 
through the vapor filter 152 (and is represented by Eq. 2 
above), is divided by the signal 155, which is the backscatter 
return from the atmosphere 144, that does not pass through 
vapor filter 152 (and is represented by Eq. 5 above). This 
removes changes in signal transmission that are independent 
of the factors to be measured: 
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pressure from the return signal. Initially, computer 156 uses 
theoretical Rayleigh functions which are lunctions of tem- 
perature and pressure in conjunction with the measured filter 
transmission to generate a lookup 170 table of laser, Ray- 
leigh, and filter (LRF(v)) convolutions that are dependent on 
atmospheric temperature and pressure. Computer 156 then 
compares the normalized return signal to the lookup table 170 
to determine atmospheric temperature and pressure. In order 
to compare the return signal with the lookup table 170 com- 
puter 156 accounts for the magnitude of the Mie scatterers as 
well as any changes in air density, which changes the magni- 
tude of the Rayleigh signal. 

The vapor filter is used as a bandstop filter for frequency 
stability, optical depth, and optimal filter shape. For the pur- 
poses of separating the Rayleigh and Mie scattering, an opti- 
cal depth of approximately 60 dB provides excellent absorp- 
tion of Mie scattering within a small frequency variance 
around v 0 (i.e., where Vyis a normalized frequency of 0 GHz). 
For example, the 60 dB of absorption in the atomic vapor cell 
exists in a region that is not contaminated by Mie scattering. 
This 60 dB absorption region is used in acquiring initial 
estimates of pressure and temperature (explained below in 
FIG. 5). Such absorption is observable below in FIG. 5 as the 
measured signal which has the magnitude of zero centered 
about v 0 . This data provides information about pure Rayleigh 
scattering that can be used to calculate the ratio of Mie scat- 
tering to Rayleigh scattering, as shown in Eq. 1 7 : 


P l T l T r E s D s J dvidv r (L(vi)F(vi - v)) 

*(v) _ (rR(v -(v r - Av D )) + mM(v -(v r - Av D ))) 
s n(v) P L T L T R E N D N fdvidv r L(vi ) 

(rR(v - ( v r - A vd)) + mM(v - (v r - Av/j))) 


(Eq 13) 


35 


fr(vo) = r E S D S yLFR(vo)+mLF(v 0 ) 
Sn( v o ) '■E^D^i r + m 


(Eq. 17) 


Since the vapor filter fully attenuates the Mie scattering in 
this region: 


Since M is a delta function, Equation 13 simplified to: 


In one embodiment, tuning laser 141 to a frequency (v^) far 
enough removed from the effects of the vapor filter 152 
enables the measurement of the ratio of the optical and detec- 
tor efficiencies of the signal channels 158 (S s as represented 50 
by Eq. 2 above) and 155 (S^as represented by Eq. 5 above). 
This enables a check for abnomalities in the filter on a scan- 
by-scan basis: 


Ss(vq) = r E S D S yLFRjvo) (Eq. 18) 

Sn(v 0 ) \-E n D n \ r + m 

where LFR(v 0 ) is the value of the theoretical return signal 
at particular atmospheric temperature and pressure. Accord- 
ingly, computer 156 calculates the ratio of Mie scattering by 
first defining a variable K 0 as follows: 


= Ss(vo) (Eq. 19) 

° S N (v 0 ) 


S s {v ) |- E S D S -| rLFR(v - Av D ) + mLF(v - Av D ) (Eq. 14) 

S N (v) VEnD N - I r + m 

45 


s s (v a ) e s d s (Eq. 15) and then solving for the ratio 

*SW( v a) E N D N 


In one embodiment, a variable can be defined as: 


m 

r 


-^-LFR{vq) - 1 


(Eq. 20) 


S s (v a ) (Eq. 16 ) 

° S N (v a ) 

Once both data sets (i.e., S s and Sjv) are symmetric about 
the same data point, computer 156 calculates temperature and 


Using the normalized signal return in the region of interest 
65 (i.e., the sloped region between the minimum and maximum 
of the signal return) and writing the result in terms of the ratio 
of m over r, yields the following: 
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S s (v) ™+^F(v) 
S N (v) a , , » 


(Eq. 21) 


5 


Substituting the ratio of m and r of Eq. 20 into Eq. 2 1 yields : 


S s (v) „ LFR(v) 


■LF(v)[l-- 


*0 


S«(») ' K{, LFR(v 0 ) +L ‘ v 'T *r„LRF(v 0 )J 
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(Eq. 22) 


Solving for LRF(v) yields: 


15 


L ^ (v) = ^™ +LF(v) [i__™ 


S N (v) K a 


K a K 0 


(Eq. 23) 


20 


where the measured signal return LRF(v) is written in 
terms of measured quantities and the theoretical values of 
LRF(v 0 ). Computer 156 then calculates LRF(v) and com- 
pares it to the lookup table 170 to determine atmospheric 2 5 
temperature and pressure, described in greater detail in FIG. 

5 . 

Accounting for power fluctuations, optical efficiencies and 
detector efficiencies, as described herein, allows for an inde- 
pendent check on vapor filter 152 during the operation OADS 30 
140 . With variable characteristics of detector channels and 
power fluctuations accounted for, computer 156 determines, 
for example, the substantially invariable characteristics of 
vapor filter 152 , such that more accurate measurements of 
received backscattered laser radiation (e.g., laser radiation 35 
148 ) are obtained. 

In one embodiment, the normalization channel depicted in 
FIG. 4 is used to remove atmospheric power fluctuations of 
laser radiation 148 . In doing so, computer 156 measures 4Q 
Rayleigh and Mie components of laser radiation 147 in terms 
of optical efficiencies and detector efficiencies. Such efficien- 
cies are typically measured on a shot-by-shot basis during the 
analysis process. In an exemplary operation, laser 141 gen- 
erates and transmits laser energy 142 as a series of pulses at a 45 
particular pulse repetition frequency (“PRF”). Computer 156 
then measures the Rayleigh and Mie components in terms of 
optical efficiencies and detector efficiencies on a pulse-by- 
pulse basis. 

To measure Rayleigh components and Mie components, in 50 
one embodiment, OADS 140 tunes the frequency of the laser 
energy 142 transmitted by laser 141 . For example, laser 141 
transmits the laser energy 142 at distal frequencies from the 
peak absorption frequency of filter 152 (illustrated by Vyin 
FIG. 4 ) to provide a frequency-independent measurement. 55 
Computer 156 then determines the line shape of laser energy 
142 through filter 152 . 

In one embodiment, measured intensity of the detected 
backscattered laser radiation (e.g., as determined by elec- 
tronic signal 158 ) is functionally compared to normalized 60 
atmospheric factors. Since the measured intensity often 
depends upon Mie scatterers (e.g., aerosols) and air density 
changes due to altitude changes and temperature changes. 
The air density changes and the temperature changes are not, 
however, removed through the normalization processes 65 
described herein; for computer 156 to accurately determine 
air parameters such as temperature and pressure of air 144 , air 
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density changes are removed from the detected backscattered 
laser radiation so that computer 156 can accurately determine 
the air parameters. 

FIG. 5 shows graph 240 with curves 241 (detected back- 
scattered laser radiation at an air density causing Rayleigh 
scattering), 242 (detected backscattered laser radiation at a 
higher air density causing both Rayleigh and Mie scattering) 
and 243 (theoretically pure Rayleigh scattering). Curves 241 , 
242 and 243 illustrate laser energy magnitudes (plotted with 
respect to axis 250 ) as a function of frequency (plotted with 
respect to axis 251 ). In one embodiment, computer 156 pro- 
cesses data from curves 241 , 242 and 243 to determine other 
air parameters. For example, Mie scattering effects are sub- 
stantially isolated and removed from calculations to deter- 
mine air temperature and air pressure, since these Mie scat- 
tering effects produce inaccurate measurements due to 
inconsistent aerosol concentrations. 

To determine the air temperature and air pressure, com- 
puter 156 processes the data from curves 241 , 242 and 243 to 
substantially isolate and remove the Mie scattering effects, 
such as those found in curve 242 . In processing the data from 
curves 241 , 242 and 243 , computer 156 calculates lookup 
table 170 in substantially real time using a measured laser/ 
filter profile (i.e., as measured at detector 165 of FIG. 2 ) 
convolved with theoretical Rayleigh functions for a particular 
temperature and pressure (e.g., illustrated by curves 242 and 
243 ). Computer 156 then scales the measured return signal 
LFR(v) (i.e., illustrated by curve 241 in this example) with the 
ratio of m to r determined by Eq. 20. Computer 156 then 
analyzes data near the deepest portion of the filter attenuation 
(i.e., approximately +/-0.5 GHz from Vy) to estimate pressure 
and/or temperature. This portion corresponds to a 60 dB 
region of absorption not contaminated by Mie scattering, a 
preferred aspect of the invention that provides temperature 
and pressure accuracy by providing a reliable temperature 
base from which to increment temperature and/or pressure 
estimates. 

Computer 156 calculates theoretical Rayleigh return 
assuming an initial temperature estimate and performs a 
Least Square Error (LSE) calculation to determine the accu- 
racy of the temperature with respect to the theoretical Ray- 
leigh function. Computer 156 repeats the process with incre- 
mental changes to temperature and/or pressure until an 
optimal fit (i.e., a LSE calculation that corresponds to design 
specifications) is achieved. Although discussed in detail with 
respect to LSE, other approximation methods, such as New- 
ton-Rhapson and Monte Carlo, may be used in alternative 
embodiments. Accordingly the invention is not intended to be 
limited to the calculation methods of the embodiments 
described herein. 

Temperature affects air density in a manner that is recip- 
rocal to pressure; increasing pressure increases density, while 
increasing temperature decreases density. Additionally, 
increasing temperature increases the Rayleigh lineshape 
width while increasing pressure increases the Rayleigh line- 
shape height. Accordingly, for each incremental value of tem- 
perature and/or pressure, the Rayleigh lineshape is unique. 
Those skilled in the art should readily understand that such 
scattering theory is generally defined in “On The Kinetic 
Model Description Of Rayleigh-Brillouin Scattering From 
Molecular Gases”, G. C. Tenti, D. Boley and R. C. Desai, 
Canadian Journal of Physics, vol. 52, pg. 285-290 (1974). 

In one example, computer 156 determines air density 
changes by aligning peak absorption frequencies of curves 
241 , 242 and 243 , illustrated at frequency \ f . Since curve 241 
represents detected backscattered laser radiation containing 
substantially no Mie scattering, curve 241 can be used as a 
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reference where Mie scatterering has been eliminated. In one 
example, computer 156, therefore, uses curve 243 to remove 
the effects of Mie scattering, aligning curves 241, 242 and 
243 and by calculating a ratio of the detected backscattered 
laser radiation to a theoretically pure Rayleigh scattering 
which then determines air density (the ratio of curves 241 and 
242). Mie scattering effects are then removed by subtracting 
curve 243 from the calculated ratio of curves 241 and 242. 
With the Mie scattering essentially removed from the mea- 
surement, computer 156 more accurately determines air tem- 
peratures and air pressures. 

FIGS. 6 and 7 show other exemplary graphs that may be 
used in determining air pressure and air temperature. FIG. 6 
illustrates a graph 260 of electronic signals 163 and 167 of 
FIG. 2 respectively generated by detectors 162 and 165 of 
FIG. 2. Graph 260 shows the electronic signals 163/167 rep- 
resenting light intensity as a function of normalized signal 
strength (axis 261) versus frequency (axis 262). FIG. 7 illus- 
trates a graph 280 of electronic signals 158 and 155 of FIG. 2 
respectively generated by detectors 153 and 154 and repre- 
senting light intensity as a function of normalized signal 
strength (axis 281), versus frequency (axis 282). These four 
light intensities (represented by electronic signals 163, 167, 
158 and 155) are measured, over time, through a collection of 
laser pulses. For example, the laser eneigy 142 of FIG. 2 
generated by laser 141 at a certain PRF is swept in transmitted 
frequency such that each laser pulse is emitted at a different 
frequency. Electronic signals 163 and 167 therefore illustrate 
how laser 141 may sweep laser energy 142 in frequency 
across an absorption band 263 of the vapor filter 152. Illus- 
tratively, FIG. 6 shows one complete frequency sweep of laser 
energy 142 generated by laser 141 and detected by detectors 
162 and 165. Similarly, electronic signals 155 and 158 of FIG. 
7 show detected signals of detectors 162 and 165 (backscat- 
tered) as laser 141 performs a complete sweep in frequency of 
laser energy 142 across the absorption band 283 of the vapor 
filter 152. 

From signals 163 and 167, computer 156 may for example 
determine a normalized filter transmission, by dividing dis- 
crete points of electronic signal 167 by corresponding dis- 
crete points of signal 163 . Similarly, computer 156 may deter- 
mine a normalized atmospheric return though vapor filter 152 
by dividing discrete points of signal 158 by corresponding 
discrete points of signal 155. These discrete points, described 
herein, correspond to individual pulses of laser energy 142. 

Using normalized calculations of filter transmission (e.g., 
from graph 260) and the normalized calculations of atmo- 
spheric return (e.g., from graph 280), computer 156 deter- 
mines relative optical efficiencies in the vapor filter 152. 

In one embodiment, computer 156 determines optical 
transmission for vapor filter 152 using the frequency inde- 
pendent components of data from graph 260; FIG. 6 (there is 
substantially no change in amplitude for signals 163 and 167 
at frequencies greater in magnitude than ±18 GHz from 0 
GHz illustrated at points 264, 265, 266 and 267). Computer 
156 therefore determines the ratio of optical transmission for 
vapor filter 152 by calculating a ratio of signal 167 to signal 
163, via frequency corresponding points of the signals, for 
points representing frequencies greater in magnitude than 
±18 GHz from 0 GHz. 

Similarly, computer 156 determines a magnitude of inten- 
sity of atmospheric -returned laser energy received through 
vapor filter 152 using the frequency independent parts of the 
data from graph 280 of FIG. 7 (there is substantially no 
change in amplitude for signals 155 and 158 at frequencies 
greater in magnitude than ±18 GHz from 0 GHz illustrated at 
points 284, 285, 286 and 287). Computer 156 thereby deter- 
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mines a ratio of atmospheric return with the laser power 
measurement by calculating a ratio of signal 158 to signal 155 
via frequency corresponding points of the signals for points 
representing the frequencies greater in magnitude than ±18 
5 GHz from 0 GHz. 

In one embodiment, computer 156 calculates a ratio of 
signal 158 to signal 155 for frequencies between ±0.5 GHz 
(illustrated at points 288 and 289). Such a frequency range 
substantially no Mie scattering of laser energy 142 for air 144 ; 
to it thus indicates substantially pure Rayleigh scattering. Com- 
puter 156 thus compares a Rayleigh to Mie scattering strength 
based upon the ratio of signal 158 to signal 155. Computer 
156 determines Rayleigh to Mie scattering strength by com- 
paring the ratio of signal 158 to signal 155 at frequencies 
1 5 between ±0.5 GHz to the ratio of signal 158 to signal 155 at 
frequencies greater than ±18 GHz from 0 GHz. In one 
embodiment, computer 156 performs similar calculations for 
“non-scattered” laser energy 142 (e.g., component 143B of 
FIG. 2) based on data illustrated in FIG. 6 using points 268 
20 and 269. Such a process is further described in FIG. 11. 

The ratios determined for the non-scattered laser energy 
142 and for the scattered laser energy 142 may be used in 
tandem to numerically calculate Rayleigh-Laser-Filter con- 
volution from data. The Rayleigh-Laser-Filter convolution is 
25 in turn compared to a look up table of theoretical Rayleigh- 
Laser-Filter convolution values to determine temperature and 
pressure. 

FIGS. 8 and 9 represent errors in pressure and temperature 
measurements obtained during experimentation. FIG. 6 
shows the percentage error in pressure measurements as a 
function of pressure and temperature, obtained using syn- 
thetic data with a 0.5% noise component. FIG. 9 shows the 
percentage error in temperature measurements as a function 
of pressure and temperature obtained using synthetic data 
with a 0.5% noise component. 

FIG. 10 shows one OADS transceiver 300. OADS trans- 
ceiver 300 is for example suitable for use as transceiver 110 of 
FIG. 1. OADS transceiver 300 is illustrated in FIG. 10 as a 
4Q midsectional view to reveal interior components of the trans- 
ceiver. In one embodiment, OADS transceiver 3 00 is hermeti- 
cally sealed and configured for projecting laser radiation 301 
to air 302 and for detecting backscattered laser radiation 303 . 
Optical coupling 305 connects to optical fiber 306. Optical 
45 fiber 306 communicates laser energy radiated from a laser to 
project laser radiation 301 to air 302. OADS transceiver 300 
may include beam steering optics 307 and lens 308 to assist in 
the projection of radiation 301. 

OADS transceiver 300 receives backscattered laser radia- 
50 tion 303 when laser radiation 301 backscatters from air 302. 
Transceiver 300 receives backscattered laser radiation 303 
through a receiving aperture 309 embedded with barrier 310. 
Aperture 309 captures and directs backscattered laser radia- 
tion 303 to rigidly mounted mirrors 311 (mounted within 
55 transceiver 3 0 0 by rigid mounting brackets 3 2 0) for improved 
detection. Mirrors 311 may be annular focusing mirrors con- 
figured for demagnifying the backscattered laser radiation 
3 03 received through aperture 3 09 . The rigid mounting brack- 
ets substantially prevent jitter of backscattered laser radiation 
60 303. 

Demagnified laser radiation 312 is reflected towards fold- 
ing mirrors 313, which, in turn, direct laser radiation 312 to 
partial reflector 314. Partial reflector 314 reflects a portion of 
laser radiation 312 in the form of laser radiation 315 towards 
65 detector 316. Partial reflector 314 is a beam splitter that 
reflects a portion (e.g., 50%) of laser radiation 312 and trans- 
mits a portion (e.g., 50%) of laser radiation 312. 



US 7,760,339 B2 


15 

Detector 316 is similar to detector 154 of FIG. 2. As such, 
detector 316 converts received laser radiation 315 into an 
electronic signal, such as electronic signal 155 of FIG. 2, for 
computer processing. A computer, such as computer 156 of 
FIG. 2 processes the electronic signal to remove laser power 
fluctuations caused by atmospheric changes in air 302. 

Laser radiation 312 that is directed to vapor filter 317 is 
filtered by vapor filter 317. In one embodiment, vapor filter 
317 is an atomic vapor filter or a mercury vapor filter and 
electronically controlled by electronics 318. Electronics 318 
may control temperature and pressure of the vapor within 
filter 317 to stabilize or electronically “lock” frequency char- 
acteristics of filter 317. Such stabilization makes OADS 
transceiver 300 less susceptible to vibrational misalignment. 

Laser radiation 312 is filtered by filter 317 and passes to 
detector 319 in the form of filtered laser radiation 321 . Detec- 
tor 319 receives and converts the filtered laser radiation 321 
into an electronic signal, which is also processed by the 
computer to determine the desired air parameters. 

FIG. 11 shows a flowchart of one exemplary methodical 
operation 400 of an OADS. A laser (e.g., laser 141 of FIG. 2) 
sweeps laser energy across a predetermined frequency spec- 
trum, in element 401. The laser may sweep the laser energy 
across a frequency range of about +/-20 GHz by firing the 
laser energy at a certain PRF. In one embodiment, the PRF is 
about 1 kHz, with a pulse width between about 50 ns and 100 
ns, and the swept frequency range is centered about a fre- 
quency corresponding to a peak absorption frequency (e.g., 
260 nm) of a vapor filter (e.g., vapor filter 152, FIG. 2). 

Laser energy is typically split into a four distinct paths such 
that the laser energy can be detected as four different inputs, 
in element 402. These four paths of laser energy correspond 
to: 1) laser energy transmitted by the laser (e.g., component 
159 of FIG. 2); 2) laser energy transmitted by the laser 
through the vapor filter (e.g., component 164 of FIG. 2); 3) 
laser energy transmitted by the laser into the air and backseat - 
tered (e.g., component 150A of FIG. 2); and 4) laser energy 
transmitted by the laser into the air and backscattered through 
the vapor filter (e.g., component 157 of FIG. 2). For simplic- 
ity, these components are hereinafter referred to as: 1) unfil- 
tered laser energy; 2) filtered laser energy; 3) unfiltered back- 
scattered laser energy; and 4) filtered backscattered laser 
energy. 

After detecting the four forms of laser energy, a computer 
(e.g., computer 156, FIG. 2), determines normalized filter 
transmission of the vapor filter, in element 403 . For example, 
the computer in one embodiment processes the unfiltered 
laser energy and filtered laser energy by dividing the magni- 
tude of the filtered laser energy by the magnitude of the 
unfiltered laser energy. The division is performed on a pulse 
by pulse basis, where divided magnitudes of the pulses have 
corresponding frequencies. 

The computer also determines, in one embodiment, nor- 
malized atmospheric return of the laser energy, in element 
404. For example, the computer processes the unfiltered 
backscattered laser energy and unfiltered backscattered laser 
energy by dividing the magnitude of the filtered backscattered 
laser energy by the magnitude of the unfiltered backscattered 
laser energy. Again, division is performed on a pulse by pulse 
basis, where divided magnitudes of the pulses have corre- 
sponding frequencies. 

Once normalized filter transmission and normalized atmo- 
spheric return of the laser energy are determined, the com- 
puter determines signal strengths for each of the filter trans- 
mission and the atmospheric return. For example, the 
computer determines the optical transmission through the 
filter by calculating a ratio of the filtered laser energy to the 
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unfiltered laser energy at particular frequency ranges, in ele- 
ments 405 and 407. The computer similarly determines the 
atmospheric return (scattering) signal strength through the 
filter by calculating a ratio of the filtered backscattered laser 
5 energy to the unfiltered laser energy at particular frequency 
ranges, in elements 406 and 408. 

The computer also determines a signal strength ratio for the 
normalized filter transmission by dividing filtered laser 
energy by unfiltered laser energy, again on a pulse by pulse 
to basis, at frequencies greater in magnitude than about +/-18 
GHz about the peak absorption frequency, in element 407. 
The computer further determines a signal strength ratio for 
the normalized filter transmission by dividing filtered laser 
energy by unfiltered laser energy on a pulse by pulse basis at 
15 frequencies between about +/-0.5 GHz, in element 405. 
These signal strength determinations correspond to fre- 
quency ranges where Mie scattering (e.g., +/-18 GHz) and 
Rayleigh scattering (e.g., +/-0.5 GHz) are most prevalent — 
useful when combined with similar signal strength determi- 
20 nations for the normalized atmospheric return. The computer 
determines a Mie scattering signal strength ratio for the nor- 
malized atmospheric return of the laser energy by dividing 
filtered backscattered laser energy by unfiltered backscat- 
tered laser energy, again on a pulse by pulse basis, at frequen- 
25 cies greater in magnitude than about +/-18 GHz about the 
peak absorption frequency, in element 408. The computer 
also determines a Rayleigh scattering signal strength ratio for 
the normalized atmospheric return of the laser energy by 
dividing filtered scattered laser energy by unfiltered backseat - 
30 tered laser energy on a pulse by pulse basis at frequencies 
between about +/-0.5 GHz in element 406. 

With signal optical transmission for the filter and signal 
strengths for both Rayleigh scattering and Mie scattering 
determined, the computer determines a Rayleigh laser filter 
35 convolution, in element 409. For example, the computer in 
one embodiment performs a convolution of the optical trans- 
mission with the Rayleigh and Mie scattering signal strengths 
corresponding to the frequency ranges for Rayleigh and Mie 
scattering of +/-0.5 GHz and +/-18 GHz, respectively. The 
40 computer then accesses a lookup table, such as lookup table 
170 of FIG. 2, having theoretical Rayleigh laser filter convo- 
lution values to determine temperature and pressure of the air, 
in element 410. 

Advantages to the embodiments described above include 
45 improved computations of various air parameters, such as air 
speed, air temperature and air pressure, substantially regard- 
less of altitude and/or Mie scattering. Other advantages 
include a system that accurately performs in a variety of 
vibrational environments. Still, other advantages include an 
50 ability to determine temperature and pressure within a par- 
ticular region of atmosphere without a prior knowledge of the 
atmosphere. Additionally, the instrumentalities as shown and 
described above advantageously require reduced on-aircraft 
system calibrations, a more robust system in vibrational envi- 
55 ronments, and system health checks. 

Since certain changes may be made in the above methods 
and systems without departing from the scope of the inven- 
tion, one intention is that all matter contained in the above 
description or shown in the accompanying drawings be inter- 
60 preted as illustrative and not in a limiting sense. By way of 
example, those skilled in the art should appreciate that the 
OADS and the OADS transceivers, as described herein, may 
be constructed, connected, arranged, and/or combined in 
other various manners without departing from the scope of 
65 the invention. Another intention includes an understanding 
that the following claims are to cover all generic and specific 
features of the invention described herein, and all statements 
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of the scope of the invention which, as a matter of language, 
might be said to fall there between. 

What is claimed is: 

1. A method for remotely sensing air outside a moving 
aircraft, comprising: 

generating laser radiation swept over a range of wave- 
lengths with a tunable laser; 

splitting the laser radiation into a projected component and 
a control component; 

projecting the projected component from the aircraft into 
the air to induce scattered radiation; 
receiving a portion of the scattered radiation as backscat- 
tered radiation; 

filtering the backscattered radiation and the control com- 
ponent with a notch filter to form a filtered backscattered 
radiation and a filtered control component; 
converting the filtered control component into one or more 
electronic control signals; 

converting the filtered backscattered radiation into one or 
more electronic backscattered radiation signals; and 
processing the one or more electronic control signals and 
the one or more electronic backscattered radiation sig- 
nals to derive a measurement of Rayleigh scattering 
from the electronic backscattered radiation signals; and 
using the measurement of Rayleigh scattering to determine 
one or more air parameters selected from the group 
consisting of an air pressure, an air density, and an air 
temperature. 

2. The method of claim 1 wherein the backscattered radia- 
tion is not passed through an interferometer prior to conver- 
sion to the one or more electronic backscattered radiation 
signals. 

3. The method of claim 1, further comprising using the 
backscattered radiation to determine a velocity of the aircraft 
through the air, and wherein neither the measurement of an air 
pressure, an air density, or an air temperature nor the deter- 
mination of velocity is performed by passing the backscat- 
tered radiation through an interferometer. 

4. The method of claim 3, wherein the backscattered radia- 
tion is measured in at least two axes and wherein an angle of 
attack is determined. 

5. The method of claim 4, wherein the backscattered radia- 
tion is measured in at least two axes and wherein an angle of 
sideslip is determined. 

6. The method of claim 1 wherein the notch filter is an 
atomic vapor filter, the electronic control signal correspond- 
ing to absorption of the filtered portion of the control com- 
ponent being used to determine absorption characteristics of 
the atomic vapor filter and the electronic backscattered signal 
corresponding to the filtered portion of the backscattered 
radiation being used to derive the measurement of Rayleigh 
scattering. 

7. The method of claim 6 further comprising: 
determining a wavelength difference between a laser wave- 
length associated with peak absorption of the control 
component in the atomic vapor filter and a laser wave- 
length associated with peak absorption of the backscat- 
tered radiation in the atomic vapor filter; and 

determining a velocity from the wavelength difference. 

8. A method for remotely sensing air parameters compris- 
ing: 

generating laser radiation, the laser radiation being tuned 
across a range of wavelengths; 
splitting the laser radiation into a projected component and 
a control component; 

projecting the projected component from an aircraft into 
the air to induce scattered radiation; 
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receiving a portion of the scattered radiation as backscat- 
tered radiation; 

passing at least a portion of the control component through 
a first path through a filter to produce a filtered control 
5 component, the filter having a peak absorption wave- 
length within the range of wavelengths; 
passing at least a portion of the backscattered radiation 
through a second path through the filter to produce fil- 
tered backscattered radiation; 

l o converting the filtered backscattered radiation into an elec - 

tronic filtered backscattered radiation signal; 
converting the filtered control component into an elec- 
tronic filtered control component signal; 
determining a frequency shift of the backscattered radia- 
15 tion by comparing a peak absorption of the electronic 
filtered backscattered radiation signal to a peak absorp- 
tion of the electronic filtered control component signal; 
and 

determining a vector component of air velocity from the 
20 frequency shift. 

9. The method for remotely sensing air parameters of claim 
8, wherein the filter is an atomic vapor filter. 

10. The method for remotely sensing air parameters of 
claim 9, wherein the filter is a mercury atomic vapor filter. 

25 11. The method for remotely sensing air parameters of 

claim 8, further comprising: 

converting a second portion of filtered backscattered radia- 
tion into an electronic backscattered radiation signal; 
using the electronic backscattered radiation signal, and the 
30 electronic filtered control component signal to deter- 

mine a relative broadening of an absorption peak of the 
filter observed with the backscattered radiation as com- 
pared to the control component; 
using the relative broadening of the absorption peak to 
35 determine an air pressure. 

12. The method for remotely sensing air parameters of 
claim 8, further comprising: 

converting a second portion of filtered backscattered radia- 
tion into an electronic backscattered radiation signal; 

40 using the electronic backscattered radiation signal, and the 
electronic filtered control signal to determine a relative 
broadening of an absorption peak of the filter observed 
with the backscattered radiation as compared to the con- 
trol component; 

45 using the relative broadening of the absorption peak to 
determine an air temperature. 

13 . A system for determining a velocity of an aircraft in air, 
comprising; 

a tunable laser for generating laser radiation; 

50 apparatus for sweeping the laser over a range of wave- 
lengths; 

a beam-splitter for dividing the laser radiation into a con- 
trol component and an emitted component; 
an atomic vapor filter having an absorption peak within the 
55 range of wavelengths; 

apparatus for directing at least a portion of the control 
component through the atomic vapor filter onto a first 
detector; 

apparatus for directing the emitted component into the air 
60 to form backscattered radiation; 

apparatus for directing at least a portion of the backscat- 
tered radiation through the atomic vapor filter onto a 
second detector; 

a computer for receiving a signal from the first detector and 
65 a signal from the second detector, for determining a 
difference between a wavelength of the laser corre- 
sponding to peak absorption of the control component 
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by the atomic vapor filter and a wavelength of the laser 
corresponding to peak absorption of the backscattered 
radiation by the atomic vapor filter, and for determining 
a velocity from the difference. 

14. The system of claim 13 wherein the atomic vapor filter 
is a mercury vapor filter. 

15. A system for sensing of air outside a moving aircraft, 
comprising: 

at least one tunable laser for generating laser energy swept 
over a range of wavelengths; 
at least one transceiver, comprising 

transmitting optics for splitting the laser energy into a 
projected component and a control component, and 
for projecting the projected component into the air; 
receiving optics for receiving scattered laser energy 
from the air; 

at least one notch filter for filtering at least one of the 
control component and the scattered laser energy, and 
detectors for converting the control component and the 
scattered laser energy into an electronic control compo- 
nent signal and an electronic scattered component sig- 
nal, the detectors not combining the scattered laser 
energy with the control component prior to generating 
the electronic scattered component signals; and 
a computer for processing the electronic signals to distin- 
guish molecular scattered radiation from aerosol scat- 
tered radiation and for determining one or more air 
parameters based on the scattered laser energy. 

16. The system of claim 15, the transmitting optics com- 
prising free space optics. 

17. The system of claim 16, the transmitting optics com- 
prising optical fibers. 

18. The system of claim 15, the receiving optics compris- 
ing one or more annular mirrors. 

19. The system of claim 15, wherein the at least one notch 
filter is an atomic vapor filter. 

20 . A method for remotely sensing air parameters compris- 
ing: 

generating laser radiation, the laser radiation being tuned 
across a range of wavelengths; 
splitting the laser radiation into a projected component and 
a control component; 

projecting the projected component from an aircraft into 
the air to induce scattered radiation; 
receiving a portion of the scattered radiation as backscat- 
tered radiation; 

passing at least a portion of the control component through 
a first path through a filter to produce a filtered control 
component, the filter having a peak absorption wave- 
length within the range of wavelengths; 
passing at least a portion of the backscattered radiation 
through a second path through the filter to produce fil- 
tered backscattered radiation; 
converting the filtered backscattered radiation into an elec- 
tronic filtered backscattered radiation signal; 
converting the filtered control component into an elec- 
tronic filtered control component signal; 
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determining a frequency shift of the backscattered radia- 
tion by comparing a peak absorption of the electronic 
filtered backscattered radiation signal to a peak absorp- 
tion of the electronic filtered control component signal; 
5 determining a vector component of air velocity from the 
frequency shift, converting a second portion of filtered 
backscattered radiation into an electronic backscattered 
radiation signal; 

using the electronic backscattered radiation signal, and the 
to electronic filtered control component signal to deter- 
mine a relative broadening of an absorption peak of the 
filter observed with the backscattered radiation as com- 
pared to the control component; and 
using the relative broadening of the spectral line to deter- 
15 mine an air pressure. 

21. The method of claim 20 wherein the filter is an atomic 
vapor filter and the absorption peak is a spectral line of the 
atomic vapor filter. 

22 . A method for remotely sensing air parameters compris- 
20 ing: 

generating laser radiation, the laser radiation being tuned 
across a range of wavelengths; 
splitting the laser radiation into a projected component and 
a control component; 

25 projecting the projected component from an aircraft into 
the air to induce scattered radiation; 
receiving a portion of the scattered radiation as backscat- 
tered radiation; 

passing at least a portion of the control component through 
30 a first path through a filter to produce a filtered control 
component, the filter having a peak absorption wave- 
length within the range of wavelengths; 
passing at least a portion of the backscattered radiation 
through a second path through the filter to produce fil- 
35 tered backscattered radiation; 

converting the filtered backscattered radiation into an elec- 
tronic filtered backscattered radiation signal; 
converting the filtered control component into an elec- 
tronic filtered control component signal; 

40 determining a frequency shift of the backscattered radia- 
tion by comparing a peak absorption of the electronic 
filtered backscattered radiation signal to a peak absorp- 
tion of the electronic filtered control component signal; 
determining a vector component of air velocity from the 
45 frequency shift; 

converting a second portion of filtered backscattered radia- 
tion into an electronic backscattered radiation signal; 
using the electronic backscattered radiation signal, and the 
electronic filtered control component signal to deter- 
50 mine a relative broadening of an absorption peak of the 
filter observed with the backscattered radiation as com- 
pared to the control component; and 
using the relative broadening of the absorption peak to 
determine an air temperature. 

55 23 . The method of claim 22 wherein the filter is an atomic 

vapor filter and the absorption peak is a spectral line of the 
atomic vapor filter. 



